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CN I In this article a series of solutions with higher baryon numbers in the chiral quark 

J>->' soliton model are reported. The chiral quark soliton model (CQSM) is a simple quark 

■ model that incorporates the basic features of QCD, e.g. the chiral symmetry and 

^ \ its breakdown accompanied by the appearance of the Goldstone bosons. It was 

00 ' shown that the baryon number one {B = 1) solution provides correct observable as 

. a nucleon including mass, electromagnetic value, spin carried by quarks, parton dis- 

^ , tributions and octet, decuplet 5C/(3) baryon spectra. The B = 2 axially symmetric 

\ soliton solution was obtained numerically by one of us (N. Sawado). For i? > 3, in a 

' series of our papers we obtained i? = 3 ~ 9, 17 minimal energy soliton solutions with 

. point-like symmetries and S = 3 ~ 5 axially symmetric saddle-point solutions, using 

\ rational map ansatz constructed for multi-baryon number skyrmions. An interesting 

' property of the solutions is that the symmetry of the background configuration is 

^) • reflected in the degeneracy of the valence quark spectra. For instance, the resultant 

I quark bound spectra are doubly degenerate for the axisymmetric solitons and are 

Q^' triply degenerate for B = 3 tetrahedrally symmetric solitons. These results confirm 

the existence of the quark shell structure. The shells consist of the four-fold degen- 
D , erate ground state and higher levels with various patterns, which are realized by the 

I interplay of two symmetries, SU{2)l x SU{2)ji symmetry for the quarks and the 

^ ■ symmetries of the chiral fields. To obtain correct physical observable, the quantum 

^ . corrections are necessary. We shall show the quantum states of the axisymmetric 

^ I solitons within the collective quantization. Upon quantization, various observable 

■ - - ' spectra of the chiral solitons are obtained. According to the Finkelstein-Rubinstein 

constraints, the quantum numbers of the solitons coincide with the physical obser- 
vations only for B = 2 and 4 while B = 3 and 5 do not. The SU{3) extension of 
the B = 2 soliton is also studied to predict various strange dibaryon states within 
this model. 



1 Introduction 



Although QCD is generally accepted as the underlying theory of the strong inter- 
action, most low- and medium-energy nuclear phenomenology may be successfully 
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described in terms of the hadronic degrees of freedom. Investigations for deuteron 
photodisintegration and deep inelastic scattering of leptons by nuclei suggest the 
necessity of including quark degrees of freedom ^ [21 El S El • For this reason it 
was suggested that nuclear theory should be reformulated to take into account the 
underlying quark theory. However, QCD is too hard to get insight into the low- and 
medium-energy nuclear phenomenology since the coupling constant become large 
at these energy scales and one can not perform perturbation as in the other gauge 
theories. It is therefore necessary to formulate low-energy effective theories for the 
strong interaction. 

The chiral quark soliton model is one of such QCD effective theories including 
quark degrees of freedom and baryons as a chiral soliton. In 1997 Diakonov et al 
predicted an exotic state with strangeness +1 within the topological soliton pic- 
ture of baryons and remarkably such resonance state has been experimentally 
discovered recently [Z]. In a naive point of view, this new exotic baryon can be 
interpreted as a five quark bound state, opening up a new paradigm of the nuclei 
as a multi-quark bound state. 

The chiral quark soliton model (CQSM) was developed in 1980's as a low-energy 
effective theory of QCD. Since it includes the Dirac sea quark contribution and 
explicit valence quark degrees of freedom, the model interpolates between the con- 
stituent quark model and the Skyrme model jHUHlEl^DEI- The CQSM is derived 
from the instanton liquid model of QCD vacuum and incorporates the nonpertur- 
bative feature of the low-energy QCD, spontaneous chiral symmetry breaking. It 
has been shown that the B = 1 solution provides correct observable as a nucleon 
including mass, electromagnetic value, spin carried by quarks, parton distributions 
and octet SU (3) baryon spectra. 

For B = 2, the stable axially symmetric soliton solution was found in Eq. |13j . 
The solution exhibits doubly degenerate bound spectra of the quark orbits in the 
background of the axially symmetric chiral field with winding number two. Upon 
quantization, various dibaryon spectra were obtained, showing that the quantum 
numbers of the ground state coincides with those of a physical deuteron . 
For B > 3, the Skyrme model shows that the soliton solutions have discrete, crystal- 
like symmetries |161 117| . From the similarity of the chiral field action between the 
Skyrme model and the CQSM one can expect that soliton solutions in the CQSM 
have same symmetries as skyrmions with the same baryon number. Since it is too 
complicated to perform a numerical computation if one imposes such discrete sym- 
metries directly on the chiral fields, Houghton, Manton and Sutcliffe thus proposed 
remarkable ansatz, rational map ansatz, for multi-skyrmions |18j . Applying this 
ansatz to the CQSM we obtained multi-chiral quark soliton solutions with point-like 
symmetries for i? = 3 ~ 9, 17 as well as saddle-point solutions for B = 5,9 [TOl OUj . 
The solutions exhibit a large degeneracy and mass gap in the valence quark orbits. 
These results confirm the existence of the quark shell structure. Esepecially the 
large degeneracy implies that our solutions with such polyhedral symmetries may 
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be the lowest-lying configuration. 

The solitons that we obtain are classical objects and therefore must be quan- 
tized to assign definite spin and isospin to them. Quantization of the solitons can be 
performed semiclassically for their rotational zero modes. Quantizing the solutions 
with discrete symmetries is, however, a formidable task in CQSM. Thus, before em- 
barking those discrete symmetries, it will be instructive to study axially symmetric 
solutions which are much simpler Besides, considering the fact that for some 
higher baryon numbers, the ground states of the skyrmions do not agree with the 
experimental observation 22^ , the possibility that axially symmetric solutions may 
provide correct ground states can not be excluded. In fact it was found in Ref. |23] 
that the axially symmmetric BPS monopoles up to chrage five have lower energies 
than those of discrete symmetries. We therefore investigate classical and quantum 
multi-soliton solutions in the CQSM with axial symmetry up to = 5. 

Since the first prediction of the H-particle in a MIT bag model calculation |24j . 
there have been many efforts to study the spectrum of dibaryonic systems includ- 
ing strangeness. We shall apply our formulation to study these six quark states 
constituting a dibaryon and make a prediction for its mass spectra. 

In Sec. Ill the formulation of the chiral quark soliton model is introduced. In 
Sec. El we obtain axially symmetric soliton solutions for i? = 2 ~ 5 and solutions 
with the polyhedral symmetries for i? = 3 ~ 9 including saddle-point configurations. 
The relation between the symmetries of the chiral fields and the degeneracy of the 
valence quark spectra is discussed in Sec. @]from a group theoretical point of view. 
In Sec. IS] we perform zero mode quantization for the obtained classical solitons. 
Imposing the Finkelstein-Rubinstein constraints on the states, the ground states of 
the axially symmetric solitons are constructed and examined if they agree with the 
experimental observation. Conclusions and discussions are in SecEl 



2 The Chiral Quark SoUton Model : General Formalism 



The CQSM is derived from the instanton liquid model of the QCD vacuum and 
incorporates the nonperturbative feature of the low-energy QCD, spontaneous chiral 
symmetry breaking. The vacuum functional is defined by Ref. |H] 



Z 



DTTVipDip^ exp 



(1) 



where the SU(2) matrix 



U + - — —U^ with [/ = exp(iT-7r//^) 



describes chiral fields, if) is quark fields and M is the constituent quark mass. /,r is 
the pion decay constant and experimentally /,r ~ 93MeV. 



N.Sawado and N.Shiiki 



4 



The B = 1 soliton solution has been studied in detail at classical and quantum 
level in Refs. [51 1^ [TUl llll I12j . To obtain solutions with i? > 1, we shall employ 
the chiral fields with winding number B in the Skyrme model as the background of 
quarks, which can be justified as follows. 

In Eq. performing the functional integral over and il)^ fields, one obtains 
the effective action 

Ses{U) = -iNcSp IniD = -iN^logdetiD, (2) 

where iD = i^—MU"'^ is the Dirac operator. The classical solutions can be obtained 
by the extremum condition of © with respect to U. For this purpose, let us consider 
the derivative expansion of the action |25[ I26[ [T2] . Up to quartic terms, we have. 



(3) 



where = U'^d^U . The baryon number B can be calculated in terms of the 
topological charge. 



j d^xti{LiLjLk) . (4) 



Suitably adjusting the coefficients C, one can identify the first two terms of Eq. Q 
with the Skyrme model action. However, the 4th order terms tend to destabilize so- 
lutions and no stable classical solution can be obtained from the above action j25U27j . 
Nevertheless, because of their similarity, it will be justified to adopt the configura- 
tions of the solutions in the Skyrme model to chiral fields in the CQSM. 

In the CQSM, the number of valence quark is associated with the baryon number 
such that the baryon number B soliton consist ol Nc x B valence quarks. If the 
quarks are strongly bound inside the soliton, their binding energy become large and 
the valence quarks can not be observed as positive energy particles [SSlUni- Thus, 
one gets the picture of the topological soliton model in the sense that the baryon 
number coincide with the winding number of the background chiral field when the 
valence quarks occupy all the levels diving into negative energy region. 

Let us rewrite the effective action in Eq. @ as 

5eff = -iNe log det{i^ - MU""' ) = -iN^ log det {idt - if (C/^^ )) (5) 

where 

H{U"'') = -ia-V + pMU'^^ . (6) 

The classical energy of the soliton can be estimated from the quark determinant 
in Eq. © [SHll^- We introduce the eigenstates of operators, idt — H{U^^) and 
i?(C/T5), such that 

H{W^^)(t>^{x) = E^<l)^{x) , (7) 
[idt -H{U^^)) = X^,n'^^,n , (8) 
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where ^^,n = e and A^^„ = —E^ + tOn- Imposing on ^^^„ the anti-periodicity 

condition, ^'^^„(a;,T) = — ^'^^^(a;, 0), reads 

cj„T = (2n + 1)^. (9) 

The determinant in Eq. © then becomes 

(2n + l)7r^ 



det{idt-H) = ]JA^,„ = + 



T 

fj,,n /i,n 



(2n + 1)2^2 



2" 



exp(i^ I^mI^) + exp(-i|i?;.|r)) (10) 



where 

n>0 



2^2 



and the product formula for the cosine function cos(z) = n^i(l — 4z2/(2n— l)^vr2) 
has been used. Inserting ([T?H) into Eq. ©, one obtains 



2 - ,111 

where is the valence quark occupation number which takes values only or 1. 
Correspondingly, the classical energy is given by 

-Ecstatic = -Eval + -E'vac (12) 

where 

Eval = ^ nfj,\Efj,\ , i?vac = ~ 2^^^ ^ '^'^^ ' 

representing the valence quark and sea quark contribution to the total energy re- 
spectively. 

The effective action Ses{U) is ultraviolet divergent and hence must be regular- 
ized. Using the proper-time regularization scheme |321) we can write 



Sim = ^Sp (e-^^^- - e- 



-N^T r ^ r -sp .-(^^+-^) - «-(^o^+-^) 



2 7„oo 27r Ji/A^ T 



(13) 
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where Dq and Hq are operators with U = 1. The total energy is then given by 



^static [C/] = E^U] + Svac[f/] " ^vac[C/ = 1] 



(14) 



where 



H — ^ exp 

i A* I 

with 



A 



1 



and -E^gj is the valence energy of the i th valence quark. A is a cutoff parameter 
evaluated by the condition that the derivative expansion of Eq. ((T^ reproduces the 
pion kinetic term with the correct coefficient, i.e., 



f 



1^ 



l/A 



dr 



2 T 



(15) 



In this model, the constituent quark mass M is the only free parameter and we take 
the value M = 400 MeV, in which the observable of the nucleon and the delta are 
well reproduced . From Eq. H15|) and by using the values of M, , we obtain 
A ~ 637 MeV. 

For B = \, one imposes a spherically symmetric ansatz [hedgehog ansatz) 



U{r) = exp{iF{r)r ■ r) = cos F(r) + ir ■ t sin F{r) , 
with the boundary condition for the profile function F{r) 

F(0) = -7r, F(oo) = 0. 



(16) 



(17) 



Substituting the ansatz ((T^ into Eq. Q) under the boundary condition ((T7j) . one 

gets 



B 



1 

2^ 



sin^ F{r) dF{r) 



Airr'^dr = — 
dr TT 



F + 



sin2F 



1, 



The one-quark hamiltonian © becomes 

H{W'') = -ia • V + /3M(cos F(r) + i'y^r • TsinF(r)) 



(18) 



(19) 



This hamiltonian does not commute with the total angular momentum J nor the 
isospin t/2 but commute with the grand spin operator K = J + t/2. H also 
commutes with the parity operator V = ^q. Hence the one-quark eigenstates are 
labeled by the K = 0, 1, 2, • • • and the parity V = ±. The three valence quarks 



N.Sawado and N.Shiiki 



7 



occupy their lowest states = 0^ and responsible for the baryon number of the 
sohton (= 1). In this context, the baryon number is not a topological origin unhke 
the Skyrme ©model. 

Field equations for the chiral fields can be obtained by demanding that the total 
energy in Eq. (|14j) be stationary with respect to variation of the profile function 
F(r), 

oFyr) 

which produces 

S{r) sin F{r) = P{r) cos F(r), (20) 

where 

S{r) = N,Y,{^AE^) + ^ig<E^)M{E^)){^l\^H{\x\ - r)|//) , (21) 
P{r) = N,Y,{nAE,.)+sign{E^)M{E^)){fi\i-f''-f''r ■ t5{\x\ - r)|/i) . (22) 

The procedure to obtain self-consistent solutions of Eq. (PH]) is that 1) solve the 
eigenequation in the hamiltonian 1)191) under an assumed initial profile function 
Fo(r), 2) use the resultant eigenfunctions and eigenvalues to calculate S{r) and 
P{r), 3) solve H2U|) to obtain a new profile function, 4) repeat 1) — 3) until the 
self-consistency is attained. 

The calculated energy of the B = 1 soliton is -Estatic = 1192 MeV with the 
constituent quark mass M = 400 MeV. 



3 The Classical Configurations 

3.1 The Axially Symmetric Configuration 

It is known that the minimal energy solution for B > 1 is not spherically symmetric. 
In the Skyrme model the configuration with i? = 2 is axially symmetric [331 1341 1351 
l36] and can be written by 

U{x) = cos F{p, z) +iT ■ nsmF{p, z), (23) 

where 

n = (sin G(p, z) cos m^ip, sin 6(p, z) sin m^ip, cos 0(p, z)) (24) 

and mw is the winding number of the pion fields. We shall use this configuration in 
the background to obtain axially symmetric chiral quark solitons. 
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The extremum conditions for the total energy 



^static [U]=0, —, r^^static P] = (25) 



SF{p,z) ' ' 5e{p,z 

yield the following equations of motion for the profile functions, 

R^{p, z) cos G(p, z) = R^{p, z) sin G(/9, z) , (26) 
S{p, z) sin F{p, z) = P{p, z) cos F{p, z) (27) 

where 

P{p, z) = R^{p, z) sin e(p, z) + R^{p, z) cos Q{p, z) . (28) 
In terms of eigenfunction (j) in Eq. (0), R^,R^ and are given by 

i?^(p, z) = R^,,{p, z) + i2L(P, , (29) 
i?^ (/>, z) = R^,, (p, z) + (p, z) , (30) 

S'(/3, Z) = S'val(/), z) + S'vac(yO, ^) (31) 

where 

R^^^{p,z) = Y2i j dip4>i{p,(p,z)i'y5{Ticosm^(p + T2smm^(p)(l)i{p,ip,z), 

i 

X (ti cos m„93 + T2 sin m^(p)(f>f^ {p, (f, z) , 

^val(/'> ^) = X] / '^'^^»(/'' ^)n5T-30i(/?, ^) , 

^vac(P> 2;) = ^J\f{Ef,) Sgn{Ef,) / dip4>^{p, if, z)i-f5T3(l)i,{p, if, z) , 

-SvaK/?,^) = X] / dip(i)i{p,ip,z)(j)i{p,ip,z) , 
i 

5'vac(p, ^;) = ^ AA(£;^) sgn(£;^) / dipcpf^ip, if, z)(j)f,{p, if, z) , 

and subscripts, vac and val, represent the vacuum and valence quark contributions 
respectively. The boundary conditions for the profile functions were constructed by 
Braaten and Carson (37) : 

F(/3, z)— >0 as p^ + z^— >oo, 

F(0,0) = -7r, G(0,z) = | . (32) 



N.Sawado and N.Shiiki 



9 



In Fig. ^ we show the spectrum of the quark orbits in the background of chiral 
fields with winding number = 2, 4, as a function of the size parameter X. The 
axially symmetric profile functions are parameterized by X as 

F{p,z) = -vr + vr + for ^/ + < X 

= otherwise, (33) 

Q{p,z) = tan-^ (p/z). (34) 

To examine the spectrum in detail, let us consider the hamiltonian defined in Eq. ® . 
For the axially symmetric chiral field in Eq. (|23j) . this hamiltonian commutes with 
the third component of the grand spin operator and the time-reversal operator 
T. These are specifically, 

1 1 

^3 = ^3 + + ^m^Ts , (35) 
T = Z7173 • iTiTsC (36) 

where L3, (T3, and T3 are respectively the third component of orbital angular momen- 
tum, spin, and isospin operator, and C is a charge conjugation operator. The parity 
operator is defined by P = 70 for odd B, and P = 70T3 for even B. Thus, the eigen- 
values of the hamiltonian can be specified by the magnitude of K-^ and the parity 
TT = ±. We have K3 = 0, ±1, ±2, ±3, • • • for odd B, and Ks = ±i, ±f , ±|, ±|, • • • 
for even B. Since the hamiltonian is invariant under time reverse, the states of +K3 
and —K3 are degenerate in energy. 

Prom Fig-^it can be seen that for m-w = 2, the bound states diving into negative 
region are doubly degenerate with K3 = ± ^ • Thus putting Nc = 3 valence quarks on 
each of the bound levels, we have the B = 2 soliton solution. Similary, for = 3, 
we obtained the spectrum of = ±1" (double degeneracy) and = O"*" states 
diving into negative-energy region which corresponds to the B = 3 soliton solution. 
For m„ = 4, the spectrum of = i^"*" and = ib| (both doubly degenerate) 
states dive into negative region, having the B = 4 soliton solution. For = 5, the 
spectrum of K3 = =b2^(double), = ibl~(double) and K3 = O"*" states dive into 
negative-energy region, having the B = 5 soliton solution. These results confirm 
that the baryon number of the soliton is identified with the number of diving levels 
occupied by Nc valence quarks. It is interesting that the degeneracy which occurs 
due to symmetry of the chiral field reduces the number of states, making large 
shell gaps. This observation indicates that degeneracy in the valence quark spectra 
contributes to make classical energies of the soliton solutions lower. In fact, our 
B = 2 solution which is considered to be the minimum energy soliton from the 
study of the B = 2 skyrmion provides the maximum degeneracy in spectra. It will 
be interesting to study minimum solutions from this point of view. 
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Table 1: The classical observable, Mass (in MeV), the mean radius of toroid and 
the root mean square radius (in fm) for i? = 1 ~ 5. 



B 






EvslI 






-Ecstatic 


{p) 


V(r2) 


1 


173 








674 


1192 




0.785 


2 


173 


173 






1166 


2204 


0.672 


0.821 


3 


173 


173 


298 




1561 


3493 


0.659 


0.854 


4 


106 


106 


232 


232 


2727 


4753 


0.971 


1.140 


5 


145 


145 


319 


319 


409 2537 


6543 


1.048 


1.225 



The baryon number density is defined by the zeroth component of the baryon 
current [5], 

(37) 



b{x) = -^{ipjo^) = 6vai(a3) + 6vac(a3) 



where 



Nr 



- 5;AA(40))sgn(^W) / d^0Wt(p,^,z),^(°)(p,^, 



A contour plot of the baryon number density for each baryon number is shown in 
Fig. 01 It can be seen that they have toroidal in shape. 
The mean radius (p) is given by 



(P) = (P)val+ (P)^ 



(38) 



where 



(P)vai = I pdpdzdipp(f>l{p, (f, z)(t>i{p, f, z) , 

I 

(P)vac = — JZ-^^^/^^^S'^^^m) / pdpdzdLpp(f)l{p,<f,z)(f>f,{p,ip,z) 



The root mean square radius is given by 



(39) 
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where 

{r^)va\ = I pdpdzdtp{p^ + 2^)(/>|(/5, (f, z)(t)i{p, ip, z) , 

m„ ^ J 

{r^)va.c = — ^Af{Ef,) sgn{Ef,) [ pdpdzd(p{p^ + z^)4>^ {p, (p, z)(j)f,{p, (f, z) . 
These values for each baryon number are shown in Tabled 

3.2 Multi- Winding Number Configurations with Polyhedral Sym- 
metries 

In the Skyrme model it is known that minimal energy configurations with B > 3 
have discrete crystal-like symmetries jl6| rather than axisymmetry. We expect that 
configurations of the CQSM inherits the same discrete symmetry as skyrmions. 
However, it is too complicated to perform a numerical computation if one imposes 
such discrete symmetries directly on the chiral fields. Thus Houghton, Manton and 
Sutcliffe proposed remarkable ansatz for the chiral fields, rational map ansatz |18j . 
According to this ansatz, the chiral fields are expressed in a rational map as 

U{r, z) = exp{iF{r)nii ■ t) , (40) 

where 

f^R = ^^-^^^(2Re[i?(z)],2Im[i?(z)], 1 - \R{z)\^) 

and R{z) is the rational map. The complex coordinate z is given by z = tan(0/2)e*'^ 
via stereographic projection. 

Rational maps are maps from CP(1) to CP(1) (equivalently, from S"^ to 5^) 
classified by winding number. In Ref. ^Hl Houghton et al. showed that B = N 
skyrmions can be well-approximated by rational maps with winding number N . 
The rational map with winding number N possesses (2A^ -|- 1) complex parameters 
whose values can be determined so as to realize minimal energies within assumed 
symmetries of the skyrmion. We shall use this ansatz for the background chiral 
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fields in the CQSM. Their exphcit forms of the map are jl8l I38j 



R5 
Re 

Rr 
Rs 

Rg 



z"^ + 2V^iz^ + 1 
- 2V3iz^ + 1 ' 

z{z^ + 3.94z^ + 3.07) 
3.07z4 - 3.94z2 + 1 ' 
z'^ + O.m 



z2(0.16z4i + 1) ' 
7/^5^6 - 7z^ - 7/V5z^ - 1 
z{z^ + 7/\/5z4 + 7z2 - 7/^5) ' 
z6 - 0.14 



Z2(0.14z6 + 1) ' 

z(-3.38- 11.19iz^ + z^) 
1 - ll.mz^ - 3.38z8 ' 
- 5) 
-5^4 + 1 ' 

5f^/3z6 - + SiVSz"^ + 1 - 1.98^2(^6 _ ^^^4 _ ^2 _^ 



9 



3iV3z^ + 9^2 - 5i^/3 - 1.98z(-i^/3z6 + ^4 + ^^^2 _ 1)) ' 



172^^ - 187^1° + 119z^ - 1 
~ z2(zi5 + 119^10 + 187^5 + 17) • (^1) 

Field equations for the chiral fields can be obtained by demanding that the total 
energy in Eq. (|14j) be stationary with respect to variation of the profile function 
F(r), 

"-E'static = , 



which produces 



5F{r) " 

S{r) sin F(r) = P(r) cos F(r), (42) 



where 

S{r) = N,Y,MEt.)+s^g^{Ef,W{E^)){fih'>5{\x\-r)\fi), (43) 
P(r) = iV, ^(n^^(ii;^) + sign{E^)Ar{E^)){fi\ij''^'nR ■ t5{\x\ - r)\fi) .(44) 
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Table 2: Mass spectra for S = 1 — 9, 17 also for some excited states B = 5* , 9* (in 
MeV) . The data for S = 2 are taken from Ref . Jl] . The ratio of the mass -Egtatic to 
B X -Eg^^d^ are compared to that of the Skyrme model [T8] . 



B 










val 












-^static 


•^^statric 

Ours 


I static 

Skyrme 


1 


173 


















674 


1192 


1.00 


1.00 


2 


173 


173 
















1166 


2204 


0.92 


0.95 


3 


210 


210 


210 














1633 


3522 


0.98 


0.96 


4 


144 


146 


146 


146 












2628 


4378 


0.92 


0.92 


5 


123 


131 


131 


139 


210 










3265 


5467 


0.92 


0.93 


6 


120 


124 


150 


150 


206 


206 








3740 


6603 


0.92 


0.92 


7 


115 


120 


120 


120 


166 


166 


166 






4554 


7478 


0.90 


0.90 


8 


97 


97 


115 


120 


139 


139 


203 


203 




5229 


8565 


0.90 


0.91 


9 


69 


101 


104 


104 


107 


166 


166 


179 


179 


6046 


9573 


0.89 


0.906 


17 


83 


95 


95 


95 


153 


156 


157 


173 


175 












177 


178 


179 


192 


194 


194 


196 


196 




10586 


18650 


0.93 


0.88 


5* 


157 


157 


157 


232 


232 










2874 


5680 


0.95 


1.00 


9* 


99 


105 


105 


121 


142 


142 


210 


210 


210 


5700 


9742 


0.91 


0.91 



The baryon density h{x) is defined by 



^(^) = l^(V'7oV') = ftvai(a^) + Afield (a;), (45) 



where 



i ^ i 

- ^sign(4o))AA(4o))<^(o)(a.)t0(o)(rr)]. (46) 

To examine the shell structure of the quarks, we evaluate the radial density for 
the ith valence quark p^*^ (r) in which the angular degrees of freedom are integrated 
out, via. 



p»(r) = J d^J sm9de 6«(r,^,^) (47) 
with the baryon number 

B = ^J drr^p^'\r). (48) 
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The profile functions for 5 = 3 — 9, 17 are plotted in Fig. [51 In Table are the 
results for the valence quark levels as well as the vacuum sea contributions. The va- 
lence quark spectra show various degenerate patterns depending on the background 
configuration. The results of the ratio of the mass -^static to B x E^^^^^ and com- 
parison to those of the Skyrme model are also in Table. El They are qualitatively in 
agreement. Using Eq. ()46|). we estimated the baryon number density (see Fig. ISJ. 
As expected the density inherits the same symmetry as the corresponding skyrmion. 

The valence quark spectra for various B are shown in Fig. [Tj These results 
strongly suggest the existence of shell structure for the valence quarks. The spectra 
show (i) four fold degeneracy of the ground state labeled by Q and various degenerate 
pattern for excited levels labeled by ^i, • " " ) (ii) a large energy gap between the 
ground state G and the first excited level ^i. Small dispersions of the spectra are 
observed in the results. In some cases they are caused by the finite size effect of the 
basis (ex. 5 = 4). Growing the size rmax and increasing the number of the basis, 
more accurate degeneracy will be attained. 

In Fig. IHl are the results of p^^\r) for B = 3 — 9, 17. The behaviour of the 
density near the origin confirms the existence of three shells {G,Ai,A2)- Q behaves 
like "S-wave" and others like "P-jL'-wave" in a hydrogen-like atom. However most 
of the densities are nearly on the same surface and very small (not zero) near the 
origin. The plateau in the density observed at the center of the nucleus [22] can 
not be attained in our solutions. Therefore one may need to employ the multi-shell 
ansatz |4nj even in the case of light nuclei. 

The solutions that we obtained here are totally classical one. To reach the physi- 
cal observations, one need to take into account the quantum corrections. Especially, 
the absolute mass tends to be much higher than the real nucleon mass due to a lack 
of the Casimir effects- the loop corrections of the order of 0{N^). For B = 1, some 
attempts for this subject have been done in the Skyrme model |4H I42| I43j and in 
the hyblid quark soliton model [1^. In both cases the calculated values properly 
reduce the large classical energy to the physical nucleon mass. For B > 1, because 
of the lower symmetry of the chiral fields, the estimations of these effects are tedious 
task. Only one work for the B = 2 chiral soliton is reported in Ref. 0^ in which 
the configuration is restricted to SO {3). Since the quantum corrections strongly 
depend on B it is difficult to extract any information of the binding energy from 
our classical solutions. 

3.3 Classical Dynamics of Slowly Moving B = 2 Solitons 

So far our concern has been restricted to static solutions, we shall now extend our 
scope to cover dynamics. Once a multi-soliton solution is obtained, one can study the 
low-energy dynamics of the multi-solitons following the geodesic approximation . 
The existence of multi-soliton solutions means that there is no net force between the 
solitons. Manton showed that the multi-BPS monopole solution is a consequence 
of the cancellation of the repulsive force by the exchange of gauge bosons and the 
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attractive force by exchange of Higgs bosons |17j. The cancehation between the 
forces can be attained when the solution saturates the Bogomol'nyi bound 48 . For 
the Skyrme model, one finds 

E > l^^i? (49) 
e 

where e is a dimensionless constant which can be fixed experimentally. It is known 
that skyrmions do not saturate this bound and have slightly higher energies. 

If the solitons which are far apart initially are given some impact parameter and 
start to move slowly towards each other, they will trace a path close to the valley of 
the potential energy without climbing the wall of the potential barrier. The valley is 
just the parameter space of the soliton solution which is called a moduli space. If the 
solitons move slowly, or adiabatically, it will be a good approximation to consider 
their dynamics only on the moduli space (truncation). Then the time evolution of 
the solitons is approximately the geodesies on the moduli space. 

In the B = 2 rational map ansatz, by imposing axisymmetry we get 

where a is a parameter of a geodesic in the moduli space with — 1 < o < 1. For o = 0, 
one recovers the minimal energy toroidal configuration. Letting the parameter a be 
time-dependent, we can examine the adiabatic time evolution of initially far-apart 
two solitons, that is, 

U{r,z;a) = U{r,z]a{t)). (51) 

The numerical solutions for various values of a are shown in Fig. |^ It is 

interesting that the famous 90° scattering is observed as in the BPS monopole !511J 
and CP{\) solitons [SJ. However, as can be seen in the energy level shown in Fig. llUl 
the solutions have higher energies for larger a, which means the moduli motion is 
not adiabatic since energy costs to separate the two solitons. This is due to the fact 
that the chiral solitons does not saturate the so called Bogomol'nyi bound. Thus 
if the two solitons are initially far apart, they will coalesce to form a bound state, 
forming a torus shape but no scattering. 



4 Symmetry and the Degeneracy of the Quarks 

The bunch of valence spectra due to the potential with discrete symmetries has been 
observed in the study of heavier nuclear systems. In Ref. the valence spectra 
are highly degenerate because the deformation of the spherically symmetric shell 
produces large shell gaps. Thus the nuclei can be considered to be more stable than 
the spherical one. As discussed in Ref. the group theory should predict the 
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level structure of pion fluctuations. However, our problem is more complicated due 
to the presence of quarks. Before discussing it in detail, let us show how the shell 
deformation is related to the degeneracy of the spectrum. 
In general, if an eigenequation given by 

Hi;^ = Ef,^^ (52) 

is invariant under a symmetric operation R ^ g, the equation transforms as 

RHi;^ = H{Ri^^) = E^tP^ . (53) 

Therefore the states {ip^,Rip^} are degenerate in energy with E^. The set of 

eigenfunctions {V'^'^^K^ = 1; • • • > ^m) belonging to a given eigenvalue E^ will provide 
the basis for an irreducible representation of the group g of the hamiltonian : 

R^\'^ = ^4'^Dlf{R). (54) 

i 

The Dirac equation 

{i-ff'd^ - M)^{x) = (55) 

is a wave equation for fermions with Lorentz-covariance, i.e. its form has to be 
invariant under a transition from one inertial system to another one. Let the co- 
ordinates of the event be x'^ for an observer A and x ^ for an observer B. Both 
coordinates are connected by the Lorentz transformation 

X = ax or {x'Y = a'^x^ . (56) 

The equation are invariant under this transformation as 

{i-f^'d^ - M)^{x) =0 4^ {ii^d^, - M)^'{x) = . (57) 

The Dirac fields are hence transformed via 

^'(x') = V'(ax) = S{a)il>{x) = S{a)i^{a-'^x') (58) 

and also 

^(x) = S-^{a)^'{x') = S-'^{a)il)'(ax) , 

ij'{x') = S{a)ij'{a-^x') ^(x) = S(a-i)^'(ax) . 

One can thus obtain the relation for the transformation operator S{a) as 



(59) 
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which reads to 

{S{a)i-f''df,S-\a) - M)i)'{x') = . (60) 
In Eq. ()56() . performing transformation to the coordinates of the system B, 

dxt' dxi^ dx'" ^^dx'" ^ ' 

one gets 

{S{a)i-f^'S-^{a)aldl - M)i;'{x') = . (62) 

Prom the covariance of the Dirac equation in Eq. ()57p. S{a) must have the fohowing 
property 

a';^r = S-Ha)YS{a). (63) 

Let us consider the transformation law for the equation including the Skyrme 
chiral fields 

{i-f^d^ - MU"'' {x))^{x) = . (64) 
If the chiral fields have a point group symmetry G such as 

Uix') = G{a)U{x)Gia)^ , {G{a) G SUi2)i) , (65) 
the Dirac equation is invariant under the Lorentz transformation 

{x'y = a'^x'^ or x' = ax (66) 

with 

accompanying a corresponding iso-rotation 

{i-ff^df, - MU^^{x))il){x) = 4^ {i-i^'d'^ - MU^"{x'))'^'{x') = , 

with 

i;'{x') = k{a)i){x) = (5(a) x G{a))i){x) . 
The operator R corresponding to this rotation is thus defined by 

iP'ix) =def mix) = k{a)^l^{a-^x) . (70) 
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It can be shown that the operator R commutes with the hamiltonian using fohowing 
commutations : 

U^^{x)i;{x) = G~^U^''{x')G ■ R-^i^'ix') 

= G~^S-'^U^^{x')ilj'{x'), 
Ri^{x') = k^{x) ^{x) = K~^RtP{x') , 

kk-^RW^{x)^{x) = U""'{x')Rtp{x') , 

.-. [R,U''']=0, (71) 

ij^-f^dkil^ix) = i-fS''aid'iip{x) 

= i5~SV'§9/V(^) 

= iS-^-f^j^Sd'ik-^iP'ix') 
= iG-^S-^j°j^Sd'ii^'{x') 

Eri7°7^afc^(x) = i7°7'^a[,V''(x') 

.-. [i2,i7V4]=0. (72) 

which reads [R,H] = 0. 

Generally speaking, if eigenequation 

H(f>k = ek(t>k (73) 

is invariant in the group G, in terms of the symmetric operation R (e G) the equation 
becomes 

RHct^k = H{Rcj)k) = ek4>k ■ (74) 

Then, (f)k and R4>k are degenerate in energy. Thus, constructing R for each symmetry 
of the hamiltonian, one should be able to deduce the degeneracy structure of the 
spectra occurring in the valence level. 

As an example, let us examine a rotational operation for the hedgehog ansatz 
given in Eq. (|16)) . The chiral fields exhibit spherical symmetry and its infinitesimal 
spatial rotation is of the form 

x"' = a''^x^ ^ x"" = x"" + e-'^x^' . (75) 

Introducing the small angle e as 

el = -el = eu e? = -ej = ea , = -e? = es , (76) 

the Dirac field transforms as 

V'(a-^x) = iPix" - e'^^x'') 

dtp 
dxP' 

ip{x'') + e-{x X V)i:{x'') 

(1 + i{x X p) ■ e)V'(x^) = exp[iZ • e]V'(x'') , (77) 



V'(^^) + 7d;(-^;:) 
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where I = x x p is orbital angular momentum operator. An infinitesimal transfor- 
mation of the Lorentz and isorotation operators are given by 

i=l+'-T-^. (78) 

Thus the transformation of the Dirac field under the Lorentz and isorotation is given 
by (see Eq. ^) 

i;'{x) = Ripix) = S X iiP{a-^x) 

= (1 + Is • e)(l + • e)(l + iix X p) ■ e)^(x) 

~ (1 + i(^5] + K + xxp)- e)V'(x) = {l + iK- e)V'(x) . (79) 

This K is a grandspin operator and a good quantum number for the hedgehog 
hamiltonian. 

The point group symmetry has a discrete operation so that it is not possible to 
perform these infinitesimal transformation analysis. Nevertheless, one can see that 
the rotation in Eq. (\7U\i produces the degenerate eigenstates. Let us show briefiy the 
derivation of the rotation operator R for B = 3 tetrahedron and the transformation 
law for the numerical basis constructed in Eq. H173|l . The B = 3 tetrahedral soliton 
is characterized by two symmetry operations JH]: -^2 x -^2 and T^. Specifically, 
Z2 X Z2 are characterized by a following Mobius transformation: 

z^-l ^ R{z)^-j^^^R'{z) (80) 



resultantly. 



- n'niz) = Y^r^(^ {2Re[R'{z)],2lm[R'{z)], 1 - \R'iz)f ) 
1 



2Re[-l/R],2lm[-l/R], 1 - \1/Rf) 



l + l/\R\^ 

= (-ni,n2, -ns) . (81) 
The transformation operator g = exp[— z^r2] ensures 

g{T-nR)g^ = i-iT2){Tini+T2n2 + T3n3){iT2) 

= -nm + n2T2 - nsTs = t ■ n'^ . (82) 

and transforms the chiral field as 



U{x') = gU{x)g^ 



(83) 
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transforms the z and R{z) as 
iz + 1 



= z' ^ R{z) 



iR+l 
-iR+l 



= R'iz), 



and hence 



-iz + 1 
ni + m2 

, , _ 2R' _l + i{R + R)- 

2R' 1 - i{R + R) - \R\'^ 



(84) 



ni — 



n'l — in'2 



713 ^ "3 



1- liJ'P 



1 + \R 



712 



= n3 + mi - 



ns — mi , 



(85) 



which yields 

n'lTi + n2r2 + 713X3 



1 1 

-(7zi + i7i2)(Ti - ir2) + -(ti'i - i7l2)(ri + iT2) + 7I3T3 
^(713 + mi)(ri - iT2) + ^(713 - mi)(Ti + zr2) + 712X3 



= nzT\ + niT2 + 712X3 . 

.-TT 1 



(86) 



The transformation operator h = exp[— z|^(ti + T2 + T3)] = (1 — i(ri + T2 + T3))/2 
ensures 

hiT-nR)^,^ = ^(1 - i(ri + r2 + r3))(ri7ii +r27i2 + r37i3)^(l + i(Ti + T2 + r3)) 

= nsTi + niT2 + 7i2r3 = r • , (87) 

and transforms the chiral field as 

U{x') = hU{x)h) . (88) 

The Lorentz transformation operator S and the operators for the chiral fields 
{y, h} corresponding to the symmetric operations {x'Y = oj^^;'^ are given by 

^2 X Z2 : 



Td-. 





1 








\ 







-1 


















1 







v 











-1 / 


/ 


1 








o\ 












1 


















1 




V 





1 










(89) 



.vr 1 



Sh = ^^P[%;^(^23 + crsi + 0-12)] 



(90) 
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The R is defined by the direct product of these rotation operators together with the 
inverse spatial rotation for the spinor such as 



Rgil){x) = [Sg X g)ijj(ag ^x) 
Rhipix) = {Sh X h)^{aj^'^x) 



(91) 
(92) 



Applying these operators to the Kahana-Ripka basis (j) = {u, v} (for detail, see 
Appendix B) we finally obtain the following transformation laws: 



RgCPKM = i-lf -^cPk-M 

Rh(t>oo = 000 • 









Rh 


it) 
















(t>a = 




>11 + 
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Rh 
























h = 
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k = 




'22 — 




(93) 
(94) 

(95) 



= «<Pio 



^c = ^(0n 






-1 








^ 








1 































1 

2 


V3 
2 











V3 
2 


1 

2 / 



</'2-2) , 4>u = (t>20 ■ 



'^2^1) 



( \ 



\<t>v / 



t>22 + (t>2-2) 



(96) 



confirming that the B = 3 tetrahedron can exhibit triply degenerate spectra. 

The numerical computation indicates that the winding number strongly couple 
the elements with different K and correlated valence spectra occur as a result. As 
can be seen from the operator K with B = 2, the degeneracy of the valence spectra 
are explained in terms of the shape deformation (symmetry) as well as the winding 
number of the chiral fields The four-fold degeneracy of the lowest states may 
be ascribed to the chiral symmetry SU{2)l x SU{2)r of the hamiltonian. The 
degenerate structure for B > 2> will be well understood if symmetric operators of 
the hamiltonian which consist of the angular momentum, spin, isospin and winding 
number, are explicitly constructed. 



5 Zero Mode Quantization 

The solitons that we obtained in the previous section are classical objects and there- 
fore must be quantized to assign definite spin and isospin to them. Quantization of 
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the solitons can be performed semiclassically for their rotational zero modes. For 
the hedgehog soUton, because of its topological structure, a rotation in isospin space 
is followed by a simultaneous spatial rotation. For the axially symmetric soliton, 
there are five rotational zero modes by rotations of iso-degrees of freedom and spatial 
rotations. 

5.1 SU{2) sector 

Let us introduce the dynamically rotated chiral fields around the classical fields Us 
U{x,t) = A{t)Us{x')A{ty, x'' = r.][B{t)]x^ (97) 

where 

r.)[B{t)] = \TT[a^B{t)a,B{t)^], (98) 

and A{t) and B{t) are time-dependent SU{2) matrices generating an iso-rotation 
and a spatial rotation respectively. By transforming the rotating frame of reference, 
the Dirac operator with Eq. 1)9 7|) can be written as 

it) = - MU^''{x,t) 

= A{t)S{t)^^[idt + i^^^'^dk' - MUf{x') + iA^ A + iS^ S]S{t)A{t)^ 

(99) 

where 

and S{t) is the rotation operator for the Dirac field and is an angle of the spatial 
rotation. Note that the gamma matrices 7^ explicitly depend on the coordinates 
and do not transform as a contravariant vector jSJj. Substituting Eq. H1UU|1 into 
Eq. (|99|). one obtains 

i~D = A{t)S{t)^-i^[idt - H{Uf) + VLA + ^B\S{t)A{t)^ (101) 

where 

VLa = iA^A = ^n'^Ta , = iS^S +{rxp)-0 = Q%Jh (102) 

with Ja = l/2eahcl^l'^ ~ ^i^ ^ ^)a- and Q.b are the angular velocity operators 
for an isorotation and for a spatial rotation respectively. Under these operations the 
efi^ective action can be written by 

Se^{U) = S,s{Us) 

- iN, Sp [log{idt - H{Uf) + VIA + ^b)] - Sp [\og{idt - H{Uf)\ 

^ S,s{Us)-\iN,Sv\ogdS + UN,Sv\ogdodo^ . (103) 
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With the proper-time regularization, it reads 



SlfiU) = SlfiUs) r - Sp[e-'''^ - e-'^odo^^] , (104) 



where 



d = iu- H{Uf) + nA + ^B, = -iuj - H{Uf) -Qa-^b, (105) 
ddt = ^ _ 2icj(0A + ^b) - [H, Qa + ^b] - i^A + ^B? , (106) 
dodo^ =uj^ + . (107) 

Assuming that the rotation of the sohton is adiabatic, we shall expand the effective 
action S^s around the classical solution Us with respect to the angular momentum 
velocity ^a and Ob up to second order [HE] 

ab 

+ /o''a4^^|(t)1^5i(t)+/o^f,0|5(t)17Ui)] (108) 
where /o,a6's are the vacuum sea contributions to the moments of inertia defined by 

I^ab = ^iVc J^/(^m,^n,A)(n|rJm)(m|r;,|n), 

n,m 
n,m 
n,m 

Io,ab = '^N,Y,fiEm,En,A){n\Ja\m){m\Jb\n) 



2 

ra,m 



with the cutoff function /{Em, En, A) 

2A e-^-M' - e-^'/^' 



f{Em,En, A) 



El-El 

, sgn(£;^)erfc(|^^|/A) - sgn(E„)erfc(|g„|/A) 

E -E ' ^ ^ 
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Similarly, for the valence quark contribution to the moments of inertia, we have 

1 ^p (val|ra|m)(m|rb|val) 



tAA 
val,afe 



tAB _ lu (val I Tg I m)(m I Jfe I val) 

m^tval 

jBA _ Ar (val I I m)(m I Tb I val) 

^val,ab-^Vc 2^ _ , 

, , -'^val 
m^tval 

jBB _^nr ST (val|-/a|?7t)(m|Jfc|val) 

^val,ab — ^-^^c • UiUj 

m^tval 

The total moments of inertia are then given by the sum of the vacuum and valence 
„„ jAA _ jAA I jAA 

Finally, the effective lagrangian is obtained as 

^ = --E'static + -^lab^^A^A + ^ab^^A^B + 2^ab^^B^B ■ (HI) 

Theoretically, these moments of inertia can be computed using the eigenstates 
of Eq. Q. However, due to the difference of the boundary conditions between the 
initial and final states of the matrix element, the moments of inertia acquire nonzero 
values with vanishing pion fields. To overcome this problem, we make the following 
replacement [SS] : 

{n\Ja\m) ^ {n\[H{U2'),Ja]\m)/{En-Em) 

= {n\[MU^\la]\m)/{En-E^) (112) 

where la = — i(r x V)^- Unless the hamiltonian explicitly depend on the coordinates, 
the numerator vanishes with vanishing pion fields. The spurious contributions to 
the moment of inertia can be removed in this way. 

From axial symmetry of the system, following relations are derived 

r. . _ n j ^ j. tAB _ jAB _ jBA _ tBA _ 
— U, I ^ J , III ~ -'22 — -'ll — -^22 — ^ ' 

tAA _ tAA tBB _ tBB tBB _ 2tAA tAB _ tBA _ tAA /-|-|oN 

-'ll — -'22 ' -'ll — -'22 ' -'33 — -'33 ' -'33 — -'33 — ''''w-'33 • V-L-Loj 

The quantization conditions for the collective coordinates, A[t) and B{t), define 
a body-fixed isospin operator K and a body-fixed angular momentum operator L 
as 

li^'^'A + li^'^'B - - tr (vlf ^) - -Ka , (114) 
+ iI'^^'b - tr (^i?^) - -La. (115) 
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Table 3: Moments of inertia (in MeV ^). 



B 




Valence 


Sea 


Total 


B 




Valence 


Sea 


Total 


2 


tAA 

tAA 
-'33 


0.00773 
0.01141 
0.00429 


0.00363 
0.00464 
0.00125 


0.01136 
0.01605 
0.00554 


4 


tAA 
-'11 

tAA 
-'33 


0.01408 
0.04272 
0.01172 


0.00959 
0.01245 
0.00074 


0.02366 
0.05517 
0.01246 


3 


tAA 
-'11 

tAA 
-'33 


0.01231 
0.02174 
0.00594 


0.00280 
0.00384 
0.00027 


0.01511 
0.02558 
0.00622 


5 


tAA 
-'ll 

tAA 
-'33 


0.02786 
0.12124 
0.01368 


0.00716 
0.01112 
0.00007 


0.03502 
0.13236 
0.01375 



These are related to the usual coordinate-fixed isospin operator la and coordinate- 
fixed angular momentum Jq operator by transformations, 

la = -Ei[Ait)]K,, Ja = -E'^[B{t)fL,. (116) 

To estimate the quantum energy corrections, let us introduce the basis functions 
of the spin and isospin operators which were inspired from the cranking method for 
nuclei 

{A,B\ll3h,Jj3l3) = ^^^2 ^^^kM)^L-m^k:^B) 

where D is the Wigner rotation matrix. Then, we find the quantized energies of the 
soliton as 

E = -^static 

where i{i + l),j(j + 1) and ^3 are eigenvalues of the Casimir operators and J^, 
and the operator K3, respectively. 

In Table 131 are the results of our calculation of moments of inertia, 1^,1 
and with B = 2 — 5. It is instructive to compare our results with the Skyrme 
model where Un = 0.0104, Vn = 0.0163 and C/33 = 0.00709 which are corre- 
spondingly our and They are qualitatively in good agreement. 

5.2 Finkelstein-Rubinstein constraints 

If a multi-skyrmion describes atomic nuclei upon quantization, it has to be quan- 
tized as a boson or as a fermion whether B is even or odd. This requirement is 
implemented in the form of Finkelstein-Rubinstein (FR) constraints [HTj . For highly 
nonlinear theory enough to possess soliton solutions a consideration of continuity 
reduces to a concept of distinct topological sector. Finkelstein and Rubinstein state 
that an equally primitive concept as continuity is between multi-valued quantized 
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systems which can possess state functions double-valued under 2tt rotation, and 
those which cannot. Both concepts are ascribed to homotopy of the map between 
the physical space and the configuration space. Indeed for the SU{2) chiral soliton 
solution to exist, the physical space must be compactified to which defines a 
topological charge characterized by an integer 

713(3^) = n. (118) 

The FR constraints arise when the space-time is suitably compactified as 

7r4(5?7(2)) = Z2 , (119) 

which takes values only —1 or +1. This allows us to quantize the solitons as either 
a fermion or a boson. 

The FR constraints for the rational map ansatz was constructed in Ref. j22j and 
Ref. jSH] and applied to predict the ground states of skyrmions up to = 22. In 
this section, we shall apply the FR constraints for the rational map ansatz directly 
to our axially symmetric multi-skyrmions and obtain their ground states. 

Following the notation in Ref. ||58j, let 5 be a rotation by a around n followed 
by an isorotation by /? around N. Then the FR constraints can be defined by 

exp(— zan • J) ex.p{—ij3N ■ = XFR{g)4' (120) 

where 

, , f 1 if contractible /-,r.-,\ 

= I -1 otherwise. ^^^^^ 

and, J and / are space-fixed spin and isospin operators respectively. ip is the wave 
function which transforms under a tensor product of rotations and isorotations. In 
particular, a closed loop is noncontractible for odd B and contractible for even B, 
which is consistent with spin statistics. Consequently, quantum numbers / and J 
are half- integers for odd B and integers for even B. 

In order to construct the ground states for a given baryon number 5, let us 
define N(L{a, P)) as a homotopy invariant for a loop L generated by rotations by 
a and isorotations by /3. Then, for the axially symmetric rational map of degree B, 
it is given by Ref. |HH] 

N{L{a,(3)) = ^{Ba-P). (122) 

It can be shown that N (mod 2) determines if the loop is contractible or not in the 
same sense as B (mod 2). Therefore, (mod 2) gives the FR constraints for each 
generator of the symmetry group of the rational map. 

The axially symmetric rational map with degree B is given by 

R{z) = . (123) 



N.Sawado and N.Shiiki 



27 



There are two symmetric generators for this rational map. One is a rotation by a 
followed by an isorotation by /3 = Ba. Substituting it into Eq. ()122|) . one obtains 
N{L{a, Ba)) = 0. The FR constraints for this loop is thus given by 

^_i^(^L3-BK,)^ = ip . (124) 

where L3, are the third component of the body- fixed angular momentum and 
the isospin operators which are related with the space-fixed operators by orthogonal 
transformations (|116j) . The other symmetry is C2 with transformation 

This corresponds to a = ^ = vr and hence iV(L(7r,7r)) = B{B - l)/2. The FR 
constraints for this loop is 



e 



-in(Li+Ki) 



^ = (-l)^(S-i)/2^. (126) 



In the following we construct the ground states consistent with the derived FR 
constraints (|124() and H126() for B = 2 — 5 with axial symmetry. 

• 5 = 2 

We find the FR constraints 

^-i^(L3-2K3)^ = tP (127) 
^.i^^Lr+K,)^ ^ (128) 

This gives the ground state as |J, L3) |/, A'3) = |1,0) |0, 0). This is in agree- 
ment with the ground state (deuteron). 

• B = 3 

We find the FR constraints 

e-i^iL3-3K3)^ = ip (129) 
^_i^(Z,,+i^,)^ = -ij; (130) 

This gives the ground state as ||)|)|^)^)- The ground state (and first "ex- 
cited" state) of B = 3 have (/, J) = (i,i) (iHe+, fH+). Then, one can not 
identify our soliton with these observed isodoublet nuclei. 



B = 4 

We find the FR constraints 

g-i7r(i3-4K3)^ = V (131) 
^-m{L^+Ki)^ ^ ^ ^ (132) 
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Table 4: B = 2, mass spectrum up to i,j < 3, ^3 < 1. 



Classification 


ihj, h) 


Parity 


Mass [Me V] 


NN(-^Si) 


(0,1,0) 


+ 


2264 


NN^So) 


(1,0,0) 


+ 


2290 


NA{^P2) 


(1,2,1) 




2399 




(1,2,0) 


+ 


2477 


iVA(35i) 


(2,1,0) 


+ 


2528 




(0,3,0) 


+ 


2576 


AA(i5o) 


(3,0,0) 


+ 


2730 


AA(5P3) 


(2,3,1) 




2762 



This gives the ground state as |0, 0) |0, 0). The ground state iHe"*" has (/, J) = 
(0,0) . Then our sohton can be identified as the "a particle". 



B = 5 

We find the FR constraints 



-iTT{L3-5K3) 



-iTT{Li+Ki)^ _ ^ ^ 



This gives the ground state as | J) |/) = ||, |) 
with the observed nuclei of i? = 5 ; (/, J) = ( 



= V' (133) 

(134) 

This is not in agreement 



1 1' 

2' 2/ 



2' 2 



) (iHe+, iLi+). 



Thus, we conclude that for even the axially symmetric solitons are possible 
candidates of the ground states of B atomic nuclei as is the case of the deuteron 
and |He while for odd B they emerge only as excited states. 

The results of the quantized energy of the axially symmetric soliton solutions 
with i? = 2,4 are shown in Table |1] and El The study of the Finkelstein- Rubinstein 
constraints indicates that the axially symmetric solution with even B has the same 
quantum number as the physically observed nuclei. Some of the states may be 
observed in experiments. Specifically, in the B = 2, we obtained the / = 0, J = 1 
{^Si : deuteron) and I = 1, J = (^5o) solutions. The energy of '^5i is lower than 
the because ifi^ > I^j^ (see Table ISJ. The order is in agreement with the 
experimental observations. For B = A, the quantum number of the ground state 
/ = 0, J = coincide with the observation. In experiment, the lowest excited state 
also has I = 0, J = 0, which unfortunately can not be explained within our scheme. 
This state can be interpreted as fHe+n bound state or fH+p resonance state rather 
than the resonance of single |He. Our formulation for the multi-soliton is based on 
the simgle bag-like picture. We expect that for some resonance states one needs 
advanced formulation including multi-fragments such as B = (3 -|- 1), (2 + 2). We 
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Table 5: -6 = 4, mass spectrum up to i < 3,j < 5, ^3 < 1. 



Classification 


{hh ks) 


Parity 


Mass [Me V] 




(0,0,0) 


+ 


4753 




(0,2,0) 


+ 


4807 


P4) 


(0,4,1) 




4808 




(1,1,0) 


+ 


4813 




(2,0,0) 


+ 


4879 


AC S3) 


(1,3,0) 


+ 


4904 




(2,2,0) 


+ 


4934.2 




(0,4,0) 


+ 


4934.3 


2A(9P4) 


(2,4,1) 




4935 


N ^A(9P5) 


(1,5,1) 




4941 




(3,1,0) 


+ 


5025 


N ^A^S,) 


(1,5,0) 


+ 


5046 


2A(954) 


(2,4,0) 


+ 


5061 


AC S3) 


(3,3,0) 


+ 


5115 


N 3A(9P5) 


(3,5,1) 




5152 


N ^ACS^) 


(3,5,0) 


+ 


5278 



observed I = 0, J = 2 with positive parity as a first excited state and this channel 
should emerge as a higher resonance (roughly 28 MeV from the ground state) in 
experiment. For odd B, the constraint of C2 in Eq. (|126j) seems to assure the validity 
of the ansatz. Indeed, it provides the ground state as I = J = 1/2 for B = 3 and as 
/ = 1/2,, J = 3/2 for B = 5, which coincide with physical observations. This seems 
to make sense since in the minimal energy configurations with discrete symmetries, 
the solutions tend to have I = J = 1/2 due to their shell-like structure. However, 
unfortunately the constraint in Eq. (|124j) forbits such states. Consequently, the 
axially symmetric solitons with odd B emerge only as excited states. The resultant 
lowest state is ^ = 3657 MeV with I = 1/2, J = 5/2 for 5 = 3, and E = 6591 MeV 
with / = 1/2, J = 7/2 for B = 5. Experimentally, no possible candidate of the state 
I = 1/2, J = 7/2 are found. 

As stated in Sec. El in the soliton approach the absolute mass always tends to 
be higher due to the lack of the Casimir effects. Therefore the total energies is 
overestimated around 0.5 ~ 1 GeV in our calculation. The one-loop corrections and 
vacuum effects should be properly subtracted in order to estimate physical mass of 
the solutions. For B = 1, the obtained mass of the nucleon the delta is E]\[ = 1260 
MeV and E/^ = 1505 MeV respectively because of the lack of Casimir effects. With 
these values, the ground states of the i3 = 2 ~ 4 exhibit bound states but the 
ground state of the B = 5 does not. For B = 2, the state / = 0, J = 3 C'^s) 
a large binding energy. The search of this resonance is interesting subject and also 
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a thorough analysis of the Casimir effects for the toroidal solitons are much desired 
for a stability argument. 

5.3 SU{3) sector 

The SU{2)l x 5C/(2)R-invariant lagrangian has a natural extension to 5'C/(3) sector 
including strange quarks. The SU{2>)l x S'C/(3)ij-invariant lagrangian in the chiral 
quark soliton model is given by 

C = iP{i^ - MW- - m)^) (135) 

where U'^^{x) = e^'^^'^a{x)\al fi, ^j^g usual Gell-Mann matrices with Aq = 

1 1. In order to estimate the effects of the symmetry breaking of SU (3) explicitly, 
we introduce the current quark mass matrix 

m = diag(m,o, rn-o, rris) = moi + Am(l — \f?>\%)l?) (136) 

where Am = ms — mo is the mass difference between the strange- and the u-,d-quark. 

Within the proper-time regularization scheme, we fix our parameters by the 
input parameters, pion mass itIt^, kion mass rrix, and the pion decay constant in 
the following formulae |59| lOOj 

'^'^^^ = /J (137) 



47r2 Jq t 



"^0^ i ^Hr)e-^''' = < (138) 







Am. NrM dr 



("^o + ^)^^y^ ^4>{r)e-^'' = ml (139) 

where M = M + nriQ. For the damping function 0(r) we introduce the two cut-off 
Ai , A2 and the parameter c via 

</>(r) = ce{T - + (1 - c)e{T - l/Kl). (140) 

Using the values M = 400 MeV , mo = 6 MeV , = 93 MeV, m^ = 138 MeV 
and Eqs. (fnTll - lfT^ we obtain the parameter set {c = 0.76, Ai = 433.8MeV, A2 = 
1512.4 MeV}. By mx = 496 MeV and Eq. ((T^ . we determine rUg = 149 MeV. 

For the extension of the three flavor soliton with B = 2, we follow the usual 
collective coordinate approach for the hedgehog ansatz with B = 1 |61l I59j . Within 
the collective coordinate approach, the extension to SU{2>) is performed by trivial 
embedding [5^ : 

C/(..*) = AW( W-') \ )AHt), (141) 
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where A{t) is the time-dependent SU{3) collective rotation matrix and A* (t) = 
^Tr{T^BTjB^) is the spatial rotation matrix. Substituting (|14H) into the quark 
determinant and transforming the rotated frame of reference, one obtains fI5] : 

iD = - MW''{x,t) - m 

^ A{t)S{t)^o{idt - H{Uf) - HsB -nA + nB)S{t)A{t)^ , (142) 



where 



= -iA^A = hl\\a , (143) 
Qb = n%{^eabclV - i{r x V),) = n%Ja , (144) 



and 



HsB = A\t)f3Am^{l - V3X8)A{t) . (145) 

S{t) is the rotation operator for the Dirac fields. In the rotating system, the quarks 
feel the induced Coliolis forces Qa,^b and Hsb- ^a,^b are the angular velocity 
operators for the right flavor rotation and the spatial rotation. The Hsb represents 
the contribution to the hamiltonian due to the SU (3) symmetry breaking. 

We assume that the rotational velocities and the mass difference are relatively 
small and the expansion in powers of 0,a,^b and Am is rapidly convergent. Ex- 
panding the quark determinant and the (valence quark) Green function in terms of 
the collective angular velocity up to the second order and the quark mass difference 
of the first order. The effective action (0) with Eq. (|142|) can be rewrite as 

S-eff = -iN.Sp log iD ^ SesiUs) + Sr + Si. (146) 

Sr, Sj mean the real and imaginary contribution to the action, which are explicitly 
written as 

Sr = -'-N,[Splog dd^ -Splogdodo^], Si = -'-N,Splogdid^)-\ (147) 

where d = idt — H{Us) — Hsb — + ^b- The expansion for the real part with 
proper-time regularization form yields 

Sr ^ -y"^Sp(/.(r)[e-'^"'^"-e-°'°'^o^"] 

jdt -\io{i - Dss) + liXo^'Xn'A - iZ^in'h + hgln%n], 

(148) 
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and 



70 = 


-Am^ '^'^Ci sgn(E'„)erfc( 

i=l,2 n 


A,, )H7o|n), 


(149) 


tAA 


i=l,2 m,n 


Xa\m){m\Xb\n) , 


(150) 


tAB 


i=l,2 m,n 


Xa\m){m\Jb\n) , 


(151) 


tBB 


1=1,2 m,n 


Ja\m) {m\Jb\n) , 


(152) 



where ci = c, C2 = 1 — c and the cutoff function f{Em, En, A) is defined in Eq. H1U9() . 
For the estimation of the imaginary part of the action, we need the following ma- 
nipulation, 



Si 



-'-N.Splog d{d))-' 
f dt f ^Trlog 



dt f'dXTrl 



W - HjUs) - HsB -^A + ^B 
w - H{Us) - HsB + nA-^B 



2tt 



-Qa + ^B 



\w - H{Us) - HsB - Xi^A - ^b). 
With replacement uj — X{^}a — ^b) and analytical continuation for uj, we obtain 

— ^A + ^B 



dt I ^ f dXTt 



w - H{Us) - HsB 



2i 



duj' {^a-^b){H{Us) + Hsb) 

op ,9 , ,TT/TT \ . TT ^2 ' ^ = 



2TT ^ u;'^ + {H{Us) + Hsb) 
Introducing the following proper-time regularization. 



(153) 



+ {^H{Us) + HsBf 



/•oo 

/ dr exp[-r(cu'2 + {H{Us) + Hsb?)] (154) 

Jl/A2 



and expand up to first order, we finally obtain the form 



Si 



dt 



and 



4^3 



m,n 



i^ffeo = Am-j=^FiEm,En,A){n\(3Xa\m){m\Jh\n) . 
^ V 3 



(155) 

(156) 
(157) 
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The moments of inertia K^j^jK^f^ are derived from imaginary part of the effective 
action thus need no regularization. The "cut-off" function F{Era, En, A) becomes 

j7fTP TP A ^^ sgn(£;^) - sgn(£;„) 

F{Em,En, A ^ CO) = _ . (158) 

Also the valence quark contributions for the moments of inertia read 

2N 

7vai = Am — (val 1 70 1 val) , (159) 



r^A A ST (val|/3Aa|re)(n|Ab|val) 



r^B A (val|/3Aa|re)(n|Jfc|val) 

n^val 

rAA _ (val|Aa|n)(n|Aj,|val) 

n^tval 

- A'. E ■ (I"!') 

n^val 
n^val 

Finally we can construct the following form for the effective lagrangian 
L = -E,t^tic-^B[Us]n\-^7{l-D88) 

and total moments of inertia acquire due to their sum, e.g. I^^f^ = /^"^ai + 
^abo- where -^static is the self-consistent classical soliton energy and DabiA) = 
^Tr(Aa^AbA^) is a SU{3) Wigner rotation matrix. 

In the evaluations of the moments of inertia, the numerical difficulties may arise. 
Due to the difference of the boundary conditions between the initial and the final 
states of the matrix elements one may obtain the spurious nonzero values of the 
moments of inertia in the absence of background pion fields. Similar problem occurs 
in the {n\Ja\m) for all values of a and {n\\a\m) for a > 4 with our harmonic 
oscillator basis. To avoid the problem, we employ the replacement of Eq. (|112() 
for {n\Ja\m) jl5j . For the matrix element of A^, similar replacement induces an 
additional spurious term proportional to the mass difference. However, within our 
perturbative treatment in the mass difference, this procedure is justified. 
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With the standard canonical quantization formula for the collective coordinates 
we obtain the following quantization prescriptions 

EM^K - I^I'^'b - K^.Dsj) , a = 1, 2, 3, 
Ra=< -Y.b{&^A-KiD^b), a = 4, 5, 6, 7, (166) 

^B, a = 8, 

and 

^^ = -E(^^''^B-^|^f^A + ^.f^8,), ^ = 1,2,3, (167) 

j 

where Ra is the right isospin generator of SU(3), and Ki represents the generator 
of the spatial rotation. 

Due to the symmetry of the soliton, only the following elements of the moments 
of inertia survive: 

tAA _ jAA tBB _ jBB 
-'11 — -'22 ) -'11 — -'22 ) 

jBB _ ^2jAA jAB _ jBA _ _^tAA 
-'33 -'33 '-'33 —-'33 — ™33 > 

tAA _ jAA _ jAA _ jAA 
-'44 — -'55 — -'66 — -'77 ' 

- K^^ 

'^11 — ^22 ! 

Kg"" = -mKis^, 

i^4r = ^^5^ = i^6r = ^7^ ■ (168) 

Thus the hamiltonian becomes 

H = Estatic + Hq + Hi, 
7 



^0 - :;jAA^^a + l^ijAA - jAA^^^i 
■^^44 a=l ^ ^11 ^44 j=i 

1 1 V"^ 9 1 / 1 1 X „9 /X 

+ -^-TbbY.^I'^-2^Taa-Taa-Jbb)^1- (169) 

^-'ll j=i ^ -'33 -'ll -'ll 

In the evaluation of the iJ, we adopt a simple perturbative treatment with the mass 
difference Am |59j . Up to first order of the Am, the Hi is written as 



1 

-7(1 - Z)88) - yaa 



Hi = -^(i-Ds8)--^{DsiRi + Ds2R2) 

'll 



T^A T^A 7 



J 33 ^44 k=A 

The hamiltonian in Eq. (|169j) is diagonalized by using following collective wave 
functions of the nonperturbative part of the hamiltonian Hq 

^?ii.,y'nn„jjM^B) = ydhi^(-l)^+^3 

xD'-%iIs,Y'NNsiA)D'*j,,.mNsiB) ■ (171) 
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With these bases, the matrix element reduces to the integral of the three Wigner 
matrices which can be easily evaluated from the SU{3) Clebsch-Gordan coeffi- 
cients j59( I63j ■ The actual computations of the Clebsch-Gordan coefficients can 
be performed by using the numerical algorithm in Ref. |64j . 



Table 6: The various moments of inertia, with nis = 149 MeV. 

Valence Vacuum Total 

l^^GeV^] 7A9 405 11.54 

I^flGeY-^] 11.19 5.46 16.65 

/^3^[GeV"^] 4.36 1.61 5.97 

I^/[GeY-^] 1.64 1.26 2.90 

Kf^ 0.285 1.38 X 10^^ 0.285 

0.297 1.60 X 10~^ 0.298 

0.255 -1.05 X 10-3 0.254 

7[MeV] 292.0 1098.8 1390.8 



Our numerical calculations were performed with the constituent mass M = 400 
MeV. For a diagonalization problem of the Dirac hamiltonian, we used the deformed 
harmonic oscillator basis [^3] which is described in detail in Appendix A. The self- 
consistent classical mass was obtained as -Egtatic = 2406 MeV which differs from the 
result of SU{2) sector. The difference arises due to the specific choice of the cutoff 
scheme in Eqs. (|137() - (|139() . The values of various moments of inertia are listed 
in Table ini The quantized states are coupled to the multiplets {10}, {27}, {35}, 
{28} corresponding to {p,q) = (0,3), (2,2), (4, 1), (6,0) respectively. In Table [71 we 
show all the mass of the dibaryon states for the multiplets. The energy levels of the 
dibaryon states belonging to each multiplets are shown in Fig. 111! 

In a pioneering work of the SU{3) collective quantization of the chiral soliton |66j . 
it is pointed out that because of the constraint = 2 which arises from the trivial 
embedding, some states in the constituent quark model are not allowed in the soliton 
solution. Hence the state we obtained is not the lowest state in the S = —2 sector 
and the configuration of H-dibaryon may have not be an axially symmetric. On the 
other hand, (/, J) = (0, 0) channel in the S = —6 sector, corresponding to di-Omega 
riQ could have a rather deeper bound. According to the data of the B = 1 hedgehog 
analysis (in Ref. |^), we expect the binding energy about ~ 200 MeV. This state 
is rather promising as a candidate of the axially symmetric dibaryon. 

In the perturbative treatment of Eqs. (|165j) . (11691) we retained only linear terms 
for the mass difference Am and used SU{3) symmetric wave functions in Eq. (|171() 
for the ground state. For the B = 1 hedgehog case, such leading order approximation 
induces discrepancy with experiment j^. But it is possible to minimize it by 
incorporating appropriate higher-order effects in the Am [2Z|- In our estimation of 
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Table 7: Absolute mass of the dibaryon (in MeV), with nis = 149 MeV. 
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(-2 2 0) 
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3597 




(-3 1 0) 


4224 




(-3 1 0) 


3904 




(-4 1 0) 


4309 




(-3 \ 0) 


3781 




(-5 i 0) 


4393 




(-4 1 0) 


3966 




(-6 0) 


4478 



moments of inertia, in order to remove the spurious contributions, we employ ad 
hoc approximation for the matrix elements and it is justified only if we confine our 
calculations up to first order. We expect an extension of our scheme to the second 
order in Am is also feasible. 

6 Concluding Remarks 

This article reports the theoretical framework and the numerical results of the chiral 
quark soliton model for higher baryon number solutions. The model was inspired 
by the instanton liquid model of the QCD vacuum and thus incorporates the ba- 
sic features of QCD, e.g. the chiral symmetry and its breakdown accompanied by 
the appearance of the Goldstone bosons. This model provides correct observable 
as a nucleon including mass, electromagnetic value, spin carried by quarks, par- 
ton distributions and octet, decuplet SU{?>) baryon spectra. For B > \, different 
topological configurations from B = 1 are needed to obtain the minimal energy solu- 
tions. For B = 2, we employed the axially symmetric configuration which produces 
the minimal energy configuration in the Skyrme model. We also investigated the 
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-B = 3, 4, 5 solitons with axial symmetry. The solution exhibits doubly degenerate 
bound spectra of the one-quark orbits. This relatively large degeneracy confirms 
that the solutions are stable local minima. For B > 2, the Skyrme model predicts 
that the solutions have only discrete symmetries. According to the prediction, we 
studied the CQSM with the chiral fields of such platonic symmetries. The discrete 
crystal-like symmetries exhibit much complicated structure and the study of such 
configurations is rather formidable task. However the analysis becomes much simpler 
when we adopt the rational map ansatz to the chiral fields since with this ansatz the 
chiral fields are separable in polar coordinates and radial coordinate, which makes 
the numerical technique developed for B = 1 applicable to find solutions with higher 
B. We showed that the baryon densities inherit the same discrete symmetries as the 
chiral fields and obtained various degenerate bound spectra of the valence quarks 
depending on the background of chiral field configurations. Evaluating the radial 
component of the baryon density, shell-like structure of the valence quark spectra 
was observed. The group theory should predict these level structures resulting from 
the symmetry of the background potential. In fact the degeneracy of the valence 
spectra are determined by the winding number of the chiral fields as well as the 
shape deformation (symmetry) of solitons. The four- fold degeneracy of the lowest 
states may be ascribed to the chiral symmetry SU{2)l x SU{2)^ of the hamiltonian. 
To get better understanding of the relation between the quark level structure and 
the winding number or the shape deformation, further analysis will be worth to be 
done in future. 

Upon quantization, we computed zero-mode rotational corrections to the classi- 
cal energy. The study of the Finkelstein-Rubinstein constraints indicates that the 
axially symmetric solution with even B has the same quantum number as the phys- 
ically observed nuclei. Some of the states may be observed in experiments. For 
example, in the B = 2, we obtained I = 0, J = 1 : deuteron) and / = 1, J = 
(^^o) solutions. The energy of ^5*1 is lower than the ^Sq. The order is in agreement 
with the experimental observations. For B = 4, the quantum number of the ground 
state / = 0, J = coincide with the observation. For odd B, the constraint of C2 in 
Eq. (|126|) seems to assure the validity of the ansatz. Indeed, it provides the ground 
state as / = J = 1/2 for B = 3 and as / = 1/2, J = 3/2 for B = 5 due to their shell- 
like structure, which coincide with physical observations. However, unfortunately 
the constraint in Eq. (|124() forbid such states. Consequently, the axially symmetric 
solitons with odd B can appear only as excited states. 

For an SU (3) extension of the model, we adopted the collective quantization 
scheme with a trivial embedding form for the chiral fields. In order to estimate the 
effects of the quark mass difference, we performed the naive perturbative method 
in terms of the mass difference. We obtained the dibaryonic spectrum coupled to 
the multiplets of {10}, {27}, {35}, {28} . The state we obtained is not the lowest 
state in the S = —2 sector. The configuration of H-dibaryon may have not be 
an axially symmetric. On the other hand, (/, J) = (0, 0) channel in the S = —6 
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sector, corresponding to di-Omega Ofi may have a rather deeper bound. This state 
is promising as a candidate of the axially symmetric dibaryon. 

In our analysis, ah obtained states seem to be deep bound states. Consideration 
of the Casimir effect is, however, necessary to determine which states are stable. For 
the B = 1 skyrmion, the Casimir energies of the rotational and the translational 
zero modes were estimated by various authors. Their predictions for the total mass 
are around -0.5 ~ -1.3 GeV gH US ESI EH ESI • ^^r the B = 2 torus, the Casimir 
energies have not been estimated yet. The thorough analysis of the Casimir effects 
is desired in order to examine the stability. 

A Numerical Basis 
A.l Kahana-Ripka basis 

The numerical method widely used in this model is based on the expansion of 
the Dirac spinor in terms of an appropriate orthogonal basis. The Kahana-Ripka 
basis 1^28: which was originally constructed for diagonalizing the hamiltonian with 
the chiral fields of i? = 1 hedgehog ansatz is a plane- wave finite basis. The basis 
is discretized by imposing an appropriate boundary condition on the radial wave 
functions at the radius rmax chosen to be sufficiently larger than the soliton size. 
The basis is then made finite by including only those states with the momentum k 
as k < kmax- The results should be, however, independent on rmax and /cmax- 

The hamiltonian with hedgehog ansatz commutes with the parity and the grand- 
spin operator given by 

K = j + T/2 = l + a-/2 + t/2, 

where j,l are respectively total angular momentum and orbital angular momentum. 
Accordingly, the angular basis can be written as 

J3T3 ™'0"3 

With this angular basis, the normalized eigenstates of the free hamiltonian in a 
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spherical box with radius rmax can be constructed as follows: 

ijK{kr)\{KK + ^)KM) 
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Nk 



—r" 



3K+1 



( ^^max ) ^ 



n-1/2 



(173) 



(174) 



and Ak = k/{Ek + M). 

The momenta are discretized by the boundary condition jK{kirraax) = 0. The 
u, V correspond to the "natural" and "unnatural" components of the basis which 
stand for parity (—1)^ and (— 1)^"*"^ respectively. 

Let us construct the trial function using the Kahana-Ripka basis to solve the 
eigenequations in Eq. Q, 



K 



lim 



E E El" 

K=0 M=-K j=l 



0') ,,0') (r 
KM,ki KMV ^ 



'^(^k\i,ii''KM 



v^llAr,e,^)]. 



(175) 



A. 2 Deformed Oscillator basis 



Let us show the numerical analysis of the eigen equations in detail. To solve the 
eigenequation of the form, 

[-ioL ■ V + /3M(cosF(p, z) + i-f^r ■ hRsinF{p, z)]<pi_,{x) = £;^(/>^(a;) , (176) 
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we introduce the deformed harmonic oscillator spinor basis which was originally 
constructed by Gambhir et al. in the relativistic mean field theory for deformed 
nuclei [HH]- The upper and lower components of the Dirac spinors are expanded 
separately by the basis as 

where ^a{x., s,t) are the eigefunctions of a deformed harmonic oscillator potential 
Vosc{p, z) = ^Mioy + ^Mioy , (178) 



and defined by 



with 



^aix,ms) = -^0lrj(/')0n.(^)e*"'^xl (179) 



= ivi:i(v^w)'"'e-^"''^'4"i(«p/>') 

ITt — 0; Ij 2, • • •, Nrmax 

n, = l,3,---,2Ar^max + l or 0,2,---,2N, 



and 



X+ = ( J ) , X- = ( J ) (180) 



depending on if the eigenvalues of the third components of the spin rris (isospin 
rrir) takes +1 or —1. The functions, and Hn^, are the associated Laguerre 
polynomials and the Hermite polynomials with the normalization constants 



These polynomials can be calculated by following recursion relations 

= nL'^ix) - (n + m)L™_i(x) (182) 
Lr'(x)=L™(x)-L-i(x) (183) 



and 



Hn+i{x) - 2xHn{x) + 2n,Hn-i{x) = (184) 
^Hn{x) = 2nHn-i{x) (185) 
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where constants Op and can be expressed by the oscillator frequencies as 



an 



MiOp 

h 



n 



(186) 



which are free parameters chosen optimally. The Nj-^^ax and N^max st-re increased 
until convergence is attained. The parity transformation rule of is given by 



^a{p, + TT, -z; s, t) = (-ir+"^$„(p, V, z; s, t) 



(187) 
(188) 



where 

has been used. The parity is + for a; + = odd, and — for a; + = even. 

There are two sets of the complete basis for each parity. One is the natural 
basis with = 0+, 1", 2+, • • • , for odd B and K[ = i"^, 1"^, • • • for even B. 
Another is the unnatural basis with = 0", 1+, 2^, • • • , for odd B and = 
\ ' i^' 1 ) ■ ■ ■ for even B The natural basis is given by 

( ^fa{0),,j,^a{0){xAs) + ^fa{l),,x'^a{l){xAs) \ 
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The unnatural basis (j)^^^ is given by replacing, a ^ /? in Eq. H189|) . By using the 
natural and unnatural basis, the eigenvalue problem in Eq. (|176|) can be reduced to 
a symmetric matrix diagonalization problem. 

Let us calculate the matrix elements of the hamiltonian below. For the kinetic 
term o; • p , we have 

(*a(0)k -PK^/J'CO)) 

= 1. f d3^4-o|(^)^^^(,),-i..o^^^)^KI(^)^^,(^),-c> 



,1 1 1 



N' , 
1 ^n' 

(*a(0)k -PN^/J'CI)) 

(wo > : w'l = WO + 1 > 0) 
(wo < : w'l = Wo + 1 < 0) 



(190) 



(191) 



In the natural basis, quantum numbers {nz,n'^) takes values (1, 2), (3, 4), • • • for 
the upper part and (1, 0), (3, 2), • • • for the lower part. In the unnatural basis, 
{nz,n'z) = (0, 1), (2, 3), • • • for the upper part and {nz,n'z) = (2, 1), (3, 2), • • • for 
the lower part. 

For the potential term /?M(cos F{p, z) + i^5T ■ fiR sin F{p, z)) , we have 
(*a(o) xi\M COS F{p,z)\ (0) xi) 

= I pdpdzM COS F{p, z)(t}l^\p)<PnA^)(t)2\p)<Pn' A^) ^ (192) 
(*a(0) xi I MiT ■ fiR sin F{p, z) (o) xi) 

= - j pdpdzM cos Q{p, z) sin F{p, z)4>\l^\ {p)<Pn. {z)^^;^^J {p)<Pn'z (z) , (193) 
{^a{o)xi\MiT -fiR sin F(p, z) | (2) Xd) 

pdpdzM sin Q{p, z) sin F{p, ^)4?l {p)<t^n-. (^)</'L?' {p)4>n'.. (z) . (194) 



Other elements can be calculated in the same manner. 

To estimate the moments of inertia, e.g., (n|Aa|?ra) let us introduce the basis for 
the strange direction The hamiltonian for the strange quark commutes with 



N.Sawado and N.Shiiki 



43 



the total angular momentum 

Js = L3 + (195) 

because of the trivial construction of the SU{3) chiral fields in Eq. 1)141(1 . As a 
result, the deformed basis has common form with the SU{2) (see Eqs. ()179() - (|189|) ). 
except in the quantum number in Eq. (|189|) as 

a(0) = {nriTiz : even,a;o = — 1/2} 
a(l) = {nr,nz : odd, oji = J3 + 1/2} 

and 

/?(0) = {nr,nz : odd, uq = J3 - 1/2} 
/3(1) = {rir, : even, wi = J3 + 1/2} . 

The unnatural basis 0^"^ is given by replacing, a ^ /3 in Eq. ()189j) . 

B Operation R to the Numerical Basis 

In this appendix, we present the construction of an operator R to the Kahana- 
Ripka basis in detail. Since it operates on the spherical harmonics, let us write the 
spherical harmonics in terms of the complex variables z and their conjugate z as 
follows : 

z = tan -e'^P , (196) 

which can be related to the usual polar coordinates by 

1 - IzP 

0036* = j—TT , exp(iip) = ± 

1 + p| 
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According to the definition of the spherical harmonics [HHl , we obtain (up to / < 3) 
1 fY 2z _ 1 /3"l-UP _ I FY 2z 



2V27rl + |zP' 2V7rl + |z|2' 2 V 27r 1 + |z|2 ' 

(198) 



4V 27rVl + |z|2; ' 2V 2^ (l + lzP)^ ' 



y^o = \xfl(sC^y-l], (199) 



1 /3- 5 2z(l - |z|2) 1 /3-5/ 2z 



2V 27r (l + |z|2)2 ' 4V 27r Vl + |z 



TT Vl + |z|2; ' 4V 27r (l + k|2)3 



8V vr v''^i + ki^^ yi + kp 

1 /y/_/l-UP\2 Al-|zl2 



1 /3-7/ /I - U|2\2 \ 2z 



1 /3- 5- 74z^(l - |z|^) 1 /5-7/ 2z 

^3-2 = T\/ — , I ,nsn , ^3-3 = 7T1 



4V 27T (l + |z|2)3 ' SV vr Vl + |z|2 

According to the (inverse) transformation of the Z2 x Z2 

z^-- = e ^Tr-9,if^7r-(p, (201) 

z 

the spherical harmonics are transformed via 

YiraiO, if) ^ Yirr^in -9,7^-^) = (-1)'-'"11_„(0, ^) . (202) 
For the (inverse) tetrahedral transformation, we have 

z^t\^, (203) 
1 + z 



N.Sawado and N.Shiiki 



45 



Thus the transformation of the spherical harmonics is given by 

^21^ ^122 + ^^21 4^2-1412-2, 

Y20 - ^y|>22 - ^l20 + ^y|^2-2 , (205) 

12- l--^y22 + ^y2r 4^2-1 + 

>2-2 ^ -^l22 + ^¥21 - ^y|5^20 , +^^2-1 - ^l2-2 ■ 

Vgk^ V^^K^ V^K^ V^5k^ V^i^^ i^^ 

o o o 4 o 00 

1/3 1 1/5 1/5 1 1/3 

- 4 V 2^33 + 2^32 + ^y-l31 - - y-y3-l - 2^3-2 " 4 ^ 2^3-3 , 

>31 ^ -^^^33 - -^^32 + g>31 + —^30 + g^3-l " ^"^3-2 " ^^^3-3 , 
I30 - -^1^33 + ^1^31 - ^13-1 + ^1^3-3 , (206) 

13- I ^ ^^33 - ^^32 - V3, + 4^y30 - \y.-. - 4^n-2 + ^^3-3 , 
>^3-2 - \^Y,, - \y,2 + ^^1^31 - ^/f n-1 + ^^3-2 - \^Y,., , 

^3-3 ~8^33 + -^^^32 — -^-^^31 + "^^^30 " ^^^^ 3-1 + "^^^3-2 — g-^3-3 • 

The operation R to the basis becomes 

i?|(/jKM) .K5;c5^f^^(5:C^,^JM'l^-3))|^r3), (207) 

J3T3 mo-3 

where \lm)' representa the transformations of the spherical harmonics in Eqs. (|202() 
or (|2()5|) - (|2()7j) . Evaluationg R for each symmetric operation, one can get the final 
results in Eqs. (|U^ - (|U7)) . 
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Figure 1: Spectra of the quark orbits of B = 2 and B = 4 with axially symmetric 
ansatz (|23j) . as a function of the sohton size parameter X |21j . 

Figure 2: Contour plot of the profile functions F{p, z),Q{p, z), of 5 = 2,3 with 
axial symmetry |21j . 



Figure 3: Contour plot of the profile functions F{p, z),Q{p, z), of i? = 4, 5 with 
axial symmetry |21j . 
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Figure 4: Contour plot of the baryon number densities b{x) [fm ^] 1)3 7|) with axial 
symmetry |21| . 



Figure 5: Self-consistent profile functions for B = 3 — 9,17 in the rational map 
ansatz calculations PUl. 



Figure 6: Surface plot of the baryon number densities b{x) (|45|1 for B = 3 — 9 and 
the excited states B = 5*,9* HO]. 



Figure 7: Valence quark spectra for S = 1 — 9, 17 pUj . 



Figure 8: Angular averaged baryon densities of ith. valence quarks p^^\r) of B = 
3 — 9, 17, with the occupation number and the eigenvalue (in MeV) PO] . 



Figure 9: Surface plot of the baryon number densities for B = 2 configurations 
dSOl) i49j. 



Figure 10: Total energies of the B = 2 configurations |49j . 



Figure 11: The dibaryon spectra jl5j . 
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